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We consider a discrete-time dynamical process on graphs, firstly introduced in connection with 
a protocol for controlling large networks of spin 1/2 quantum mechanical particles [Phys. Rev. 
Lett. 99, 100501 (2007)]. A description is as follows: each vertex of an initially selected set has 
a packet of information (the same for every element of the set), which will be distributed among 
vertices of the graph; a vertex v can pass its packet to an adjacent vertex w only if w is its only 
neighbour without the information. By mean of examples, we describe some general properties, 
mainly concerning homeomorphism, and redundant edges. We prove that the cardinality of the 
smallest sets propagating the information in all vertices of a balanced m-ary tree of depth k is exactly 
(m*^^ + ( — l)'')/(m + 1). For binary trees, this number is related to alternating sign matrices. 



I. INTRODUCTION 

Background. In view of applications like quantum 
RAM or charge-coupled devices, Burgarth and Giovan- 
netti f3| introduced a protocol for arbitrarily control net- 
works of coupled spin 1/2 quantum particles (for exam- 
ple, an array of trapped ions). An important feature 
of the protocol lies on the ability of transforming the 
physical state of the entire network, by acting sequen- 
tially with the same local operation on a specific subset 
of particles. This is valuable, since physical operations 
on quantum objects are generally difficult to implement. 
It has been shown in [3| that a network can be prepared 
in an arbitrary state by acting on the particles of a sub- 
set, only if that subset satisfies certain conditions related 
to the eigensystem of the Hamiltonian. Such conditions 
can be lifted from the physical scenario and analyzed in 
a purely combinatorial setting. This can be described 
in what follows as a discrete-time dynamical process on 
graphs. 

Definition. Let G ~ (V, E) be a simple undirected 
graph with V{G) ^ {1, 2, n}. Given a set 5 C V{G), 
let N[S] = {w e V{G)\S -.^v S : {v,w} (z E{G)} be 
the (closed) neighborhood of S. Let Vs ■ [n] — > ViG) 
be a map associating a subset of V{G) to each time t £ 
[n] — {0, 1, n — 1}. We consider the following process: 

• We select a set S" C V{G) and fix -Ps(O) = S. 

• For each t e N\{0}, we have Vs[t) = Vsit - 1) U 
T, where T C N[Vs{t - 1)]. Moreover, if -u; e 
'Ps{t)\Vs{t - 1) then there \s v ^ Vsit - I) such 
that {v, w} e E{G) and AfH\{w} = Vs{t - 1). 

In words, at time i = 0, we select a subset of vertices 
7^5(0). At time i = 1, we may insert some vertices into 
^5(0) and obtain 7^5(1). The propagation will go on 
until eventually Vs{k) — V{G), for some k. Clearly, 
k < n — 1. However, the propagation is not free, but it 
obeys some rules. Specifically, a vertex w can be inserted 
at time t in Vs{t), only if it is adjacent to v & Ps{t — I) 
and all other neighbors of v (except w itself) are already 



in Vs{t - 1). If there is k such that Vs{k) = V{G), 
we say that S propagates to G. We denote this fact by 
Sr\G. Notice that Vsit - 1) C Vsit). We denote by 
fl^A the cardinality of a generic set A. The figure below 
represents the steps of the described dynamical process 
taking place on a small graph. The square vertices are 
the elements of Vs{Q)- The grey vertices are the elements 
of V{G)\Ps{t). In the first row of the figure, Vs{0) is on 
the left and 'Ps'(l) on the right; in the second row, Vs{2) 
andPs(3) = F(G). 
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Interpretation. We may depict the above scenario in 
a more concrete language: each vertex of an initially se- 
lected set has a packet of information (the same for ev- 
ery vertex in the set), which has to be diffused among 
the vertices of the graph; a vertex v can pass its packet 
to an adjacent vertex w only if w is its only neighbour 
still without the information. In this way, v does not 
need to discriminate among its neighbours, even if it is 
permitted to pass the information to only one of those. 
Equivalently, this can be interpreted as a procedure for 
coloring (with the same color) the vertices of a graph, 
in such a way that a vertex can be colored at a certain 
time step, only if it is the unique uncolored neighbour of 
an already colored vertex. Notice that the propagation 
is not synchronized, that is, we do not require that at a 
certain time k, a vertex is necessarily included in Vs{k) 
if it is the unique uncolored neighbor of a colored vertex. 

A quantitative question. Here is a precise mathemat- 
ical problem: given a simple undirected graph G, find a 
set of minimum cardinality that does propagate in G. 
The cardinality of such a set will be denoted by tt{G). 
Obviously 7r(G) is invariant under isomorphism, since it 
does not depend on the labeling of the vertices. Notice 
that in this problem we do not take into consideration the 
propagation time. The problem can be in fact modified 
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by imposing time constraints. Finding 7r(G) is of prac- 
tical importance, when trying to optimize the number 
of local operations required to initialize, and then con- 
trol, a networks of spin 1/2 particles. The computational 
complexity aspects of the question, a formulation as an 
orient at ionproblem, and approximation algorithms are 
studied in |]|. Roughly speaking, it looks like that tt{G) 
depends on the expansion properties of G. Intuitively 
higher is the number of "ways out" from each subset of 
vertices of a certain size and higher is tt{G). 

Structure of the paper. We will focus mainly on trees. 
Apart from the introduction, the paper contains two sec- 
tions. In Section 2, we underline some general properties 
of propagation, by taking as example paths, combs and 
stars. We focus on homeomorphism and the maximum 
possible number of edges that a graph can have, given 
the cardinality of the initial set ^^(O). We finally make 
a comment about hamiltonicity and propagation in di- 
graphs. In Section 3, we will focus on balanced trees. 
The main technical tool is a simple proof that the min- 
imum cardinality of 'Ps(O), such that S propagates in a 
balanced binary tree, realizes the Jacobsthal sequence. 
Our reference on the theory of graphs is the book by 
Diestel H. 

II. GENERAL FACTS BY EXAMPLE 

Homeomorphism. It is worth keeping in mind that 
7r(G) is invariant under homeomorphism. It is then more 
appropriate to think about 7r(G) not as a quantity as- 
sociated to a single graph G, but rather to a family of 
graphs, whose members are all the graphs homeomorphic 
to G. Recall that graphs G and H are homeomorphic if 
H (G) can be obtained by subdivision and smoothing on 
G (H): a subdivision of an edge {u,v} consists of delet- 
ing {u,v}, adding a vertex w, plus the edges {u,w} and 
{w,v}] a smoothing is the reverse operation, and it is 
then performed only on vertices of degree two. A tree is 
a graph in which any two vertices are connected by ex- 
actly one path. This property, plus the homeomorphism 
remark, make propagation on trees particularly amenable 
to quick observations. 

Paths. Let P„ be the path of length n. The path P„ 
models the classic spin chain with equal nearest neigh- 
bor couplings, or more complex networks via the notion 
of graph covering and equitable partitions (see, e.g., [^). 
This is probably the graph for which is simplest to de- 
termine 7r(G). Indeed n{P„) — 1, for every n. This 
fact is self-evident and it does not need a proof. If 
V{Pn) = {l,2,...,n}and{u,w} e E{G) only if u; = v+1, 
then it is sufficient to take 5 = {1} or 5 = {n}. Clearly, 
rs{n - 1) = V{Pn). Let G„ = {V,E) be a graph in 
which #V^(G) = n. Let Q — {G„ : G„ satisfies a prop- 
erty P for all n} be a family of graphs. We can look 
at the number 7r(G„) as a function / : Z — >Z, defined 



as f{n) = 7r(G„), for all G„ G G- We have seen that 
7r(P„) = 1 independently of n. This suggests the follow- 
ing structural graph theory problem: characterize classes 
of graphs Q, for which 7r(G„) = c, where c is a con- 
stant, for all G„ G Q. Paths have this behaviour, since 
T^{Pn) = 1, for all n. The same can be said for n-cycles, 
since 7r(G„) — 2. For complete graphs this is a linear 
function: 7r(i^„) = n — 1. The next figure illustrates how 
{l}r\P4 (in the obvious way): 

■ — • — o — o ■ — • — • — • 

Adding edges (I). Since 7r(P„) = 1 and the mini- 
mum degree of P„ is 1, it is also natural to ask about 
graphs for which 7r(G) is exactly the minimum degree, 
that is, the trivial lower bound. Paths suggest also 
another question: given a graph G on n vertices and 
S C V{G), what is the maximum number of edge that 
G can have such that Sr\G. One can obtain a cy- 
cle Cn from a path P„ by adding an extra edge to 
P„. For covering G„ by propagation, we need #5' > 2. 
Specifically, if S contains just two adjacent vertices then 
SrxCn- Can we augment G„ by extra edges and keep 
#5 = 2? Let £;(G„) = {l,2},{l,n},{2,3},...,{n-l,Ti}. 
Let S = {l,n}. StiU Sr\Cn + {2,n}. More generally, 
S'r^Gra -I- IJj {2, i}. It is plausible to conjecture that when 
#5 = 2 and V{G) = n, then #E{G) = 2n - 3 (n > 2) is 
the maximum possible number of edges that G can have 
if Sr\G. The graph G„ -I- Uii^i*} attains the bound. 
The graph G5 + {2, 5} U {3, 5} is drawn below. 




f^:: i y . i y ■ 

Combs. A comb Pn^k is a path P„ having a copy of Pk 
attached to each vertex. Usually the plane embedding 
of this tree is such that the copies of Pk are all drawn 
in the upper region of the plane determined by P„. This 
justifies the term "comb" ; the path Pn is then called bone 
and the paths Pk are called fingers. So, ^V{Pn,k) = kn. 
The comb Pn,k has 2 vertices of degree 2, n — 2 vertices 
of degree 3, and kn — n + 2 vertices and k vertices of 
degree 1. Given the invariance under homeomorphism, 
it is sufficient to deal with P„.2. In fact, longer fingers 
attached to the vertices of the bone P„ would not modify 
Tr{Pn,k)- Equivalently, 7r(P„,fe) = 7r(P„_2) for every k. We 
have 7r(P„ 2) = n/2 if n is even and 7r(P„ 2) = \n/2~\ if n 
is odd. The figure shows how a 3-element set propagates 
in P5,2- 

inniiiii 
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Adding edges (II) . As we have already seen in the pre- 
vious paragraphs, in some situations one can add edges, 
and a pre-selected set will still propagate in the graph. 
Certainly, one can always add edges connecting only the 
pre-selected vertices and create a clique of size ^S. In 
combs, differently from the case of C„, we can construct 
a clique of size n/2, containing then — j edges. This 
implies that the total number of edges is going to increase 
with a faster peace than in C„ augmented by redundant 
edges. In a generic graph, we can always add redun- 
dant edges connecting vertices in S and then complete 
the subgraph induced by 5' to a clique, without altering 
the dynamics. Other edges can be added provided that 
these satisfy some conditions. Given v G Vsit), if we can 
add {v, w] e E(G) then: 

1. Suppose w € Vs{t)- Then there is a vertex y such 
that y r\ w a,t time t' > t. 

2. Suppose w ^ Vsit). 

(a) If N[v] C Vs{t) then we can simply proceed 
to include w in Vsit + 1). 

(b) If N[v] ^ Vsit) then there is a vertex y such 
that y r\ w a.t time t' < t" , where t" is the 
time step at which v r\ z, for some vertex z. 

We ask: what is the maximum number of edges in 
Pn,2 + H so that 7r(P„,2) = 7r(P„,2 + H)7 What about 
graphs in general? The answer is not immediate and we 
leave it as an open problem. 



The figure shows propagation on the comb P3.2, saturated 
with an additional number of edges. The extra edges are 
represented by dotted lines. Note that adding a single 
vertex to a graph, in which we have already fixed the 
initially selected vertices, may be sufficient to stop the 
propagation. 

Stars. The complete bipartite graph Ki^n-i is also said 
to be a star on n vertices. Contextually to quantum 
networks, properties of free bosons hopping on star net- 
works where investigated in Q. If S* ^ Ki^n-i then 
= n — 2, by taking n — 2 leaves (i.e., the vertices of 
degree 1). Among all graphs on n vertices, the complete 
graph Kn is the only graph for which the ratio n/^VsiO) 
is smaller. If we include the root in S (i.e., the vertex of 
degree n — 1), then = n — 1. This implies that we can 
add — + 1 redundant edges to Ki^i^i and obtain 
Kn such that S r\ Ki^n-i and S r\ Kn, for exactly the 
same set S. If a graph G has Ki n~i as a spanning sub- 
graph then #Ps(0) £ {n - 2, n - 1}, for G. Recall that 
a spanning subgraph is a subgraph that contains all the 
vertices of the original graph. Valuable remarking that if 



a set propagates in a graph then it will propagates in all 
of its spanning subgraphs. Equally, the minimum car- 
dinality of such a set is nonincreasing when restricting 
ourselves to spanning subgraphs. 

Digraphs and hamiltonicity. An orientation of a 
graph G is a directed graph G obtained by giving a di- 
rection to the edges of G and in this way substituting 
E{G) with a set of directed arcs. The propagation dy- 
namics induces a partial ordering on the vertices of G 
and therefore an orientation. Observe that the defini- 
tion given in the introduction of this paper can be ex- 
tended to digraphs in a straightforward way. We ask: 
given a graph G, can we always find an orientation G 
and a set S such that S r\ G and #5 = 1. If G 
has a Hamilton path then the answer is in the affirma- 
tive, because we can just orient forward the edges of the 
Hamilton path and backwards the remaining edges. No- 
tice that the grid considered in 0] is Hamiltonian. We 
can then always take a single vertex to propagate in an 
arbitrary large grid, as far as we give a proper orienta- 
tion to the edges. Formally, let P„ be the Hamilton 
directed path and let 5 = {!}. Within respect to G, we 
have Vsit) ~ {l,2,...,t}, for each t. The arc set of G 
is E(G) = E(P„) U {(u, v) : {u, v} £ E{G) Av < u}. 
In the oriented version, we get #Ps(0) — 1 also for the 
complete graph. For Kn, the orientation giving rise to a 
Hamilton directed path is obtained by constructing any 
Hamiltonian tournament. Here is a picture of a single 
vertex propagating in an orientation of K^: 




III. BALANCED TREES 

Let us denote by T2,k a balanced binary tree of depth k. 
Then#V^(r2,fc) = 2*^-1. The root of r2,fc is denoted by w. 
All the remaining vertices are denoted by v^, where x S 
{0, 1}', for i = 1, k — \. In particular, {v, vq\, {v, vi\ € 
E{Tk,n) and {vx,Vxo}, {vx,Vxi} e E{T2,k), for every x € 
{0,1}* and i = l,...,k — 2. The number of time steps 
needed to cover a tree is necessarily equal to the diameter 
of the tree. 

Top-down propagation. A set S is said to propa- 
gate in T2^k by topdown propagation, when it propa- 
gates in T2,k and if Vx G 'Ps(t) then x G {0, 1}\ with 
i < t + 1, for every t. Equivalently, in topdown propa- 
gation, the information flow goes from the root to the 
leafs. As a consequence, v G Vs{0). It is straight- 
forward to determine T^TD{T2,k), ^-S-, the cardinality of 
the smallest set covering T2,fc by topdown propagation: 
given a tree T2^k, we have TTTD{T2,k) = 2^^"^. Let us 
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see why. If S* = {v} then ■Ps(l) = S. So, we need to 
include vq or vi in S. Let us take S — {v^vq}. Now, 
Psil) = {w,fo,wi}- However, Vsi2) = T'sll)- So, we 
need to include iioo, I'd, I'lo or vu in 5. In this way, if 
S = {v, vo}U{v^ : X = yOAy e {0, 1}\ 1 < i < k - 2} 
then Sn.G. It is then clear that ■KTD{T2,k) = #5* = 2'=-i. 

Bottom-up propagation. A set S is said to propa- 
gate in T2.k by bottom-up propagation, when it propa- 
gates in T2,fc and if Vx £ "Psit) then x G {0, 1}*, with 
z = A; — lift = 0, i = — 2ift = l, and so on. Equiv- 
alently, in bottom-up propagation, the information flow 
goes from the leafs {i.e., the vertices of the form v^, with 
X G {0, 1}'^"^) to the root. It is obvious that T2,k is cov- 
ered by bottom- up propagation if 7^5(0) is the set of all 
leaves, that is #Ps(0) = 2^"^. This is not equal to the 
optimum, if we want that v E S only if w is a leaf, with- 
out any further constraint. The figure shows bottom-up 
propagation in T2.4. 





We will give a formal proof of the next result, a technical 
lemma for establishing Theorem [21 

Lemma 1 For T2.k be a balanced binary tree of depth k. 
Then 

2^ + (-1)''""^ 

T^leaji^^.k) — T^iT2,k) — ^ , 

the (fc — l)-th Jacobsthal number. 

Proof. First, take T2.3. It is useful to write 7^1(0) and 
when considering T2^k- Suppose the elements of ■Pl(O) 
being leaves only. We can start by including in a single 
vertex, say vqq, in agreement with the notation defined. 
We will think of the information flow going from bottom- 
left to bottom-right, with the propagation starts from vqq 
and ending at wn. (Just think of T2^3 drawn on the plane 
in the usual way.) We will add vertices in S online, as 
required, every time the propagation stops. In this way, 
we provide that #5* is as small as possible. We have 

Woo ^ "^0 only if vqi G "PKO). The notation is easy: wqo 
propagates to vq at time 1. Now — {voo,voi}. So, 

2 1 3 A 1 • 4 1 .r 

Vq 1^ V and V r\ vi. At this stage, vi r\ vn only it 
?;io G S^. At the end = {woo, ^'oi, I'lo} and #5*^ = 3. 



It follows that TTieaf{T2,3) = 7r(T2,3) = 3. The tree 
T2^4 can be constructed by taking two copies of T2^3 and 
adding an extra vertex adjacent of the roots of the smaller 
trees. The new vertex is the root of T24. We can include 
the set in ^1(0) for T2,4. The root of T2^4 is auto- 
matically covered and so the vertex vi G V{T2^4). This is 
sufflcient to show that the set S"^ = {vooo, I'ooi, '''loO) viio} 
propagates in T2a- So, #P|(0) = (2 • #P|(0)) -1 = 5. 
Let G be the graph obtained by adding a pendant ver- 
tex to the root of T2A {i.e., a vertex of degree one). For 
this graph tt{G) — #7^1(0) + 1 = 6 puts in evidence 
a recursive way to obtain TTieaf{T2,k)- Since r2.A;+i is 
constructed with two copies of T2^k plus a new root, we 
can write #7^1(0) = 2(#7'|"^(0)) + 1, for k odd, and, 
#7?|(0) = 2{#Vs~\0)) - 1, for /c even. Such quantities 
are exactly Trieaf{T2,k) — 7r(T2,fc), because the vertices 
inserted online in S'^ are leaves. Odd and even cases 



are combined in the formula 7r(r2,fc) = 
{k — l)-th Jacobsthal number [8i] 



_ 2'' + (-l) 
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It has been pointed out that a set rv G if rep- 
resents a configuration incompatible with a nontrivial 
eigenstate of the network Hamiltonian 0] . Connections 
between the ground state vector for some special spin 
systems and the alternating-sign matrices (ASMs) form 
an active field of research in the interface between combi- 
natorics, statistical mechanics and condensed matter (see 
[3| and the references contained therein). The number of 
ASMs of size n is A{n) — YI^JIq '^^n+iy ■ ^^'^Y and Sellers 
[3] proved that A{n) is odd if and only if n is a Jacobsthal 
number. This observation could reveal a potential link 
between ASMs and the physics {e.g., properties of the 
eigensystem) involved in the protocol proposed in [3] , for 
networks modeled by trees. 

Theorem[2]is the main point of this section. The proof 
is essentially the same as the one of Lemma [TJ The se- 
quences realized by 7r(T,„_fc) can be seen as generaliza- 
tions of the Jacobsthal numbers. 

Theorem 2 LetT„i^k be an balanced m-ary tree of depth 
k. Then 



'!^leaf{Tm,k) = T^{Tm,k) = 



rn'^+i + (-1)" 



m + 1 

The table below contains the first values of n{Tm,k)'- 



m/k 
1 
2 
3 
4 
5 
6 



1 1 



3 
1 
3 
7 
13 



4 

5 

20 
51 



21 104 



5 
1 

11 
61 
205 
521 



6 


21 
182 
819 
2604 



31 185 nil 6665 



The number 



fc+i 



Qra,k=Y.^-lf^'-'m^+^- 



i=0 
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is the mfc-th entry of the table, disregarding of the signs. 
All sequences reahzed by the rows of the table appear 
to count walks of length k between any two vertices 
in the complete graph Km+i, *-e-, these are equal to 
A{Kjn+i)ij (with i ^ j), where A{Km) is the adjacency 
matrix of Km- It is an open problem to exhibit a bi- 
jection between each element in an initially selected set 
of minimal cardinality propagating in Tm,k and walks of 
length k in Km+i- 
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